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Abstract 

In the context of the Minimal Supersymmetric extension of the Standard Model, 
we analyze the production at future high-energy e^e~ colliders of second and third 
generation scalar leptons as well as scalar quarks in association with neutralinos 
and charginos, e^e~ — > f fx- In the case of third generation squarks, we also 
discuss the associated production with gluinos. We show that the cross sections for 
some of these three-body final state processes could be significant enough to allow 
for the detection of scalar fermions with masses above the kinematical two-body 
threshold, ^/s = 2m j. We then discuss, taking as a reference example the case of 
scalar muons, the production cross sections in various approximations and make a 
comparison with the full four-body production process, e^e~ f fxXi i^i particular 
around the two-sfermion threshold. 
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1. Introduction 



The search for Supersymmetric (SUSY) particles is one of the major goals of future high- 
energy colliders. Because of their strong interactions, scalar quarks and gluinos can be 
best searched for at hadron machines such as the Tevatron and the LHC 0, where they 
are copiously produced. On the other hand, the weakly interacting particles, the left- and 
right-handed scalar leptons, and as well as the charginos and the neutralinos, Xi,2 
and Xi-45 can be more efficiently probed at high-energy e^e~ colliders where the signals 
are cleaner and the backgrounds are smaller 0, ||. The cleaner environment also allows 
for the detailed study of the properties of these particles and provides the possibility to 
reconstruct parts of the SUSY Lagrangian at the low energy scale, opening a window for 
the determination of the structure of the theory at the very high energy scale 0. 

In the Minimal Supersymmetric extension of the Standard Model (MSSM) |]^, |^, the 
lightest SUSY particle (LSP) is the lightest neutralino Xi- If R-parity is conserved, this 
particle is stable and since it is electrically neutral, it is invisible and escapes experimental 
detection [at least in the simplest pair production processes]. In models where the gaugino 
masses are unified at the GUT scale, as in the minimal Supergravity model (mSUGRA) 
[^, the masses of the lightest charginos and neutralinos are such that: m^o ~ m^± ~ 2m^o 
in the case where they are gaugino-like or m^o ~ m^± ~ m^o in the case where they are 
higgsino-like. Thus, the states Xo, and Xi are not much heavier than the LSP and might 
be the first SUSY particles to be discovered. In models with a common scalar mass at 
the GUT scale, sleptons (in particular ^R) are in general much lighter than squarks and 
can have masses comparable to those of the lighter neutralinos and charginos. In the case 
of the f slepton, mixing effects that are proportional to the mass of the partner fermion, 
generate a possibly large splitting between the mass eigenstates, making that one of them 
is much lighter than the other sleptons. [The mixing in the case of selectrons and smuons, 
as well as for first and second generation squarks, is extremely small and can be safely 
neglected.] The scalar quarks of the third generation, Q = t and 6, because of the large 
Yukawa couplings of their partners quarks, can also mix strongly leading to one mass 
eigenstate which is rather light. These particles can therefore be among the first ones to 
be accessible at the next generation of e^e~ colliders. 

In e^e~ collisions, left-and right-handed scalar leptons can be produced in pairs 

e+e- ii£* {i = L, R) (1) 

if the center of mass energy is high enough, ^/s > 2m^. In the case of second and third 
generation sleptons, the production mechanisms proceed through s-channel gauge boson 
exchange, 7 and Z-boson for charged sleptons and only Z-boson for the sneutrinos. Mixed 
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pairs of left- and right-handed charged sleptons [there is no right-handed sneutrino in the 
MSSM] cannot be produced in e+e~ colhsions, since 7 and Z bosons do not couple to 
states. For first generation sleptons, additional channels are provided by the t-channel 
exchange of neutralinos in the case of selectrons and charginos in the case of sneutrinos, 
making the situation slightly more complicated. For third generation squarks, the pair 
production of the lightest states, e^e~ —>■ titl and bibl, proceeds only through s-channel 
7 and Z-boson exchange |TDI . 



Far above the kinematical threshold, the production of charged sfermions [ignoring the 
t-channel exchange diagrams] is dominated by the photon exchange and the cross sections, 
normalized to the QED point like cross section for muon pair production ctq = 47ra^/(3s), 
are simply given by Rj ~ N^Q'^^/A, where is the electric charge and Nj the color factor 
(= 3 for squarks and 1 for sleptons) of the sfermion /, and the factor 1/4 comes from the 
absence of spin sum for scalar particles compared to the case of final fermions. The cross 
sections are therefore rather large: in the case of smuons, one has a{e'^e~ f^tlM) ~ 100 
fb at a cm. energy i/i = 500 GeV, leading to a sample of the order of 50 000 events 
with the high integrated luminosities, / £ ~ 500 fb~^, expected at future linear colliders 
[as is the case for the TESLA machine [Q for instance] . This large number of events will 
therefore allow very detailed studies of the properties of these particles. 

Due the P-wave nature of the production mechanisms for scalar particles [we again 
ignore the t-channel contributions], the cross sections are strongly suppressed by = 
(1 — 4mjr/s)^/^ factors near the kinematical thresholds. Scanning near these thresholds 
is very important since it allows us to determine accurately the masses of the produced 
sfermions [^. In this case, a very refined analysis of the cross sections ||TT], |12[, including 
the non-zero decay width of the produced states and some higher order effects [such as 
Coulomb re-scattering effects and initial or final state radiation], has to be performed. 

An interesting question to ask is: what if the center of mass energy is slightly below 
the kinematical threshold for sfermion pair production, y/s < 2m j? In this case, one has 
to consider the production of off-shell sfermions, subsequently decaying into neutralinos 
or charginos and fermions. In fact, in the general case, one has to consider the associated 
production of one sfermion, its partner fermion and a light neutralino or a chargino 
state [and only in the case where they are lighter than the sfermion as is the case of the 
neutralino LSP]. These are three-body production processes and the cross sections should 
be suppressed by an extra power of the electroweak coupling oew and by the virtuality of 
the exchanged (s)particles. However, since the pair production cross sections are rather 
large, one might hope that the three-body cross sections are not too small if one is not 
far above the two-body kinematical production threshold. 
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These higher order processes are worth studying for at least four good reasons: 

(i) They provide the possibihty to discover sfermions even if they are too heavy to be 
produced in pairs, i.e. they increase the discovery reach of the e+e~ colhders. 

(ii) In the case where the sfermion is the next-to-hghtest sparticle, it will decay 100% of 
the time into its partner fermion and the LSP and the information on the sfermion- 
fermion-neutralino coupling is lost; reducing the energy below threshold where one 
has only the three-body process, allows an access to this coupling. 

{in) To make a threshold scan for the measurement of the masses one has to include the 
non-zero decay widths of the sfermions, i.e. to consider the four-body production 
process of two fermions and two gauginos; it is then interesting to compare the cross 
section in this more complicated case with the one of the three-body production 
process, i.e. with only one resonant sfermion, which is easier to handle. 

(iv) In the strong interaction sector, it provides access to the gluino which, if it is heavier 
than squarks, cannot be produced otherwise in e+e~ collisions [except via loops and 
the cross sections in this case are too small]. 

In this paper, we therefore analyze, in the context of the MSSM, the associated pro- 
duction of sfermions with neutralinos and charginos (and in the case of third generation 
squarks, with gluinos) at future high-energy e^e~ colliders. For sleptons, we consider the 
case of left- and right-handed smuons, stau's (for which the mixing will be included), as 
well as their partners sneutrinos; the more complicated case of associated production of 



e and Ue with gauginos will be discussed in a forthcoming paper []13|. For squarks, we 
will discuss the associated production of the lightest states ti and 61 with neutrahnos, 
charginos as well as with gluinos. The amplitudes for these processes will be given and 
the cross sections for masses above the two-body kinematical threshold will be shown. 
In the case of fin, flL and z/^, we will also discuss in details some approximations and the 
behavior of the cross sections near the kinematical thresholds. 

The remainder of the paper is organized as follows. In the next section, we discuss the 
sfermion-gaugino associated production mechanisms and present the amplitudes for the 
various channels. In section 3, we give some numerical illustrations for the production of 
second and third generation neutral and charged sleptons as well as for third generation 
squarks, including the associated production with gluinos). In section 4, we discuss various 
approximations and the threshold behavior of the cross sections. A short conclusion is 
given in section 5. In the Appendix, we present the analytical expression of the differential 
cross section for the associated production in the very good approximation where only 
the s-channel contributions are taken into account. 
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2. Production mechanisms and amplitudes 

2.1 The Feynman diagrams for the various processes 

The Feynman diagrams contributing to the associated production in e^e~ collisions of 
squarks and second/third generation sleptons with the lighter chargino xf ^^nd neutralinos 
X?2) that we will call electroweak gaugino^ for short, are displayed in Fig.l. More 
explicitly, the various possibilities which will be studied in this paper are, for sleptons 

fi'fl^'^Xi,2 > l^f^i^Xi with i = L,R 
^^^^Xi,2 > ^t'^rXi with i = 1 

i^rr'^Xi , i^Ti^TXi,2 (2) 
and for squarks [the diagrams for gluino final states are shown in Fig. 2] 

e+e" ^ iiixi,2 > ^i^Xl , hig 

e+e" hibx% > ~^itxt , ^ihg (3) 

In the case of smuon final states, the production is mediated by s-channel 7 and Z 
boson exchanges with either a smuon (la) or a muon (lb) is virtual and goes into a 
(charged/neutral) lepton or slepton and a (neutral/charg)ino. Additional contributions 
come from the production of chargino and neutralino pairs, either through s-channel 
gauge boson exchange (7, Z for xf and only Z for Fig. Ic) or through t-channel first 
generation slepton exchange: Ve for xf , Fig. Id or Id' (depending on the isospin of the final 
sfermion) and cl^r for final states. Fig. Id and Fig. le (note that here, there are both 
t- and M-channels because of the Majorana nature of the neutralinos) with the virtual 
gaugino going into a smuon/neutrino or smuon/muon final states, respectively. Note that 
in Fig. Ic, there is no additional diagram where the slepton-lepton pair is emitted from 
the other neutralino since the gauge boson-x°-x° vertex is already symmetrized. 

In the case of sneutrino final states, the same diagrams in Fig. 1 contribute except for 

the fact that there is no 7 exchange in diagrams (la) and (lb). For fi final states, an 

additional complication comes from the mixing which has to be taken into account in the 

xf -f-i/,- Xi'^^"'^ ci^id Z-f-f vertices and the exchange of the heavier f2 slepton in Fig. la 

when the Z boson is exchanged. 

^In the range of parameter space where the production cross sections are favorable, the hghter charginos 
and neutrahnos wiU be indeed gaugino-hke as wiU be discussed later. We will not consider the case of the 
associated production with the heavier chargino and neutralinos since the cross sections will be suppressed, 
at least by the smaller phase-space. Furthermore, for third generation sleptons and for squarks, we will 
consider only the more favorable cases of the lighter fi as well as ii and 61 states. 
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The case of associated production of third generation squarks with their partners 
quarks and electroweak gauginos is to be handled similarly as the case of f sleptons, i.e. 
to include the mixing in the various couplings and the exchange of the heavier squark 
states in diagram (la) when the Z-boson is involved. For the associated production 
of squarks and gluinos, the situation is relatively simpler since only two diagrams are 
involved: squark and quark pair production with the virtual squark or quark going into 
a pair of quark/gluino or squark/gluino; Fig. 2. 





Figure 1: Generic Feynman diagrams contributing to the associated production of 
sfermion-fermion-gaugino final states in e+e" collisions. 




Figure 2: Feynman diagrams contributing to the associated production of squark-quark- 
gluino final states in e'^e" collisions. 
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In the previous processes, all neutralino and chargino final states have to be included, 
provided their masses are such than m^^. < y/s — mj — nif, i.e. when the phase space 
is large enough for the associated production to take place. For the exchanged gauginos 
in Figs. Ic-le, all states should in principle be included. However, if at the same time 

y/s > rriy.. +my.^ and m^^ > mj+nif, both the two-body production process e+e~ XiXj 
and the two-body decay Xi ~^ ff can occur, and the total cross section would be then 
simply the cross section for the two-body production process times the branching ratio 
for the decay channel. Since this situation has been already analyzed in several places, 
we will not discussed it further here. 

The channels contributing to the various processes discussed above are summarized in 
Table 1, following the labels of Fig. 1. In the case of the associated production of sfermions 
with charginos, there is an overall minus signal difference between the contributions of 
the the s-channel gauge boson exchanges (Fig. Ic) in the case of isospin +1/2 and —1/2 
because of the charge of the chargino. In addition, diagram Fig. Id' contributes only to 
dhilhjuii'ijxt filial states, while diagram Fig. Id contributes to both UL{yi)d(l)Xi and 
//x° final states; with the conventions adopted, the common structure of the matrix 
elements for the latter two processes differs only by a relative sign. 



Final states 
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Table 1: Contributing diagrams to the various final states. 



The relative signs among different diagrams, when contributing to the same final 
state, arising out of the anticommuting nature of the fermionic fields (Wick's theorem) 
are summarized in Table 2 below. 



Diagrams 




b^ 




az 


bz 


Cz 


d 


d' 


e 


Relative Signs 


+ 


+ 


+ 


+ 


+ 


+ 




+ 





Table 2: Relative signs among different contributing diagrams. 
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Finally, we note that the production of the charged conjugate states has also to be 
taken into account. Due to CP-invariance, these cross sections are the same as for the 
corresponding previous ones, which have thus to be multiplied just by a factor of two. 



2.2 The amplitudes of the various contributions 

The generic process considered here, including the momenta pj's of all particles, is 

e+(pi) e~{p2) fj{p3) Xi{P4) f'{P5) (4) 

where x stands for an electroweak gaugino X12 xf- The amplitudes of the various 
channels, following the labels of Fig. 1, are given in a covariant gauge by: 

= e'Q/ [^e(Pl) ih- A- h)Ue{p2)] [%i (^4) ^/(Ps)] 



M2 



Sw „ (Pi + P2Y {(P4 + Pbf - m) } 

Jn 

e^Qf [VejPl) T Ue{p2)] [UxAPa) G^" {^3+ ^4 + mf) 7a Vf{p5)] 

Sw (Pi + P2y {{P3 + P4)^ - rnl} 



" ~ Sw ^ (Pl+P2)M(P3+P5)^-m2±} 

j^Z ^ ^^ZUh y [^e(Pl) (i>3- A" A) (c^Pl + C^rPr) {^^XPa) ^/(Ps)] 

sWw n {(Pl+P2)'-m|}{(p4+P5)'-"^H 

= ^ gV(^yv^) [vMriclPL + c^nPR)ueiP2)] 

{{Pi + P2Y - m^z) l(P3 + PaY - mj} 

xKM Gtr {h+ h + mf) 7a (c{Pl + c^Pr) Vf{p,)] 



W-w j {{Pi + P2Y - ^z} {iP3 + Pr.Y - ml.} 

la gZxiXj (i^3+ A ^Xj) GcZ '^fiPb)] 
^ ^ [^e(pi) (A+ A - "^X.) Gf'^'^Vfips)] [u^M Gdlln) Ue{p2)] 



et^^ {(P3 +P5)2 - ml.} {{P2-PaY - Ke(.,H),i>e)) 

^ _ [^e(pi) ggj V(P4)] [%(P5) G^^^" (i'3+ A + m^) u^] 

~ ^Wj {{Ps+PsY -ml±} {{Pi-P^Y -TTT-IJ 

e3 [v^{p,) G^l^^) V(P4)][%(P5) Gt^^-{h+ h + m^^)Gti^^R) Ue{p2)] 

s'wl^/r {iP^+P^Y-ml.}{{p,-pAY-m\,J 



The left- and right-handed projectors are Pl,r = |(1 =F 75)- The constants and the 
couphngs used in writing down these amphtudes are described and defined below, in 
terms of the electromagnetic coupling constant = Ana, the sine and the cosine of the 
Weinberg angle, sw = sm9w and cw = cos9w, T^f^{f) and Q f^(^f)i which respectively, 
are the third component of the weak isospin and the charge of the fc-th chiral sfermion 
(fermion). Note that in the case of gluino final states, only a subset of diagrams will 
contribute and the QCD factors have to be included. 

Fermion and sfermion couplings with gauge bosons: 



a 



w 

zuf, - dkMikMjk with ak = T^j^ - Q si, (6) 

k=L,R 

The 2x2 matrix M is the orthogonal rotation matrix connecting the mass and the 
chiral (L, R) eigenstates of the sfermions and defined as 

a \ / J \ 

(7) 

where, 6j is the mixing angle. For the first two generations of sfermions the mixing effect 
can reasonably be neglected (i.e. 9j = 0) so that M/„ = if Z 7^ n. Particularly, for 
sneutrinos, only Mil 7^ 0. 
Gaugino-gaugino-Z boson couplings: 

We followed Figs. 75(b,c,d) and eqs. (C88) of Haber and Kane 0] for the 'jxtxj (with 
Q ±,0 = ±1,0 being the charge of the on-shell gaugino) and the ZxiXj couplings, the 
latter being defined as gzxiXj ~ ^ij^L + O^jPr, with: 






For charginos 



1 1 

For neutralinos : = -Og = --Ni^Nj^ + -N^Nji (8) 



Fermion- sfermion- gaugino couplings: 

For the couplings of charginos and neutralinos to fermions and scalar fermions we fol- 



low Figs. 22,23 and 24 of Haber and Gunion |T^. In the amplitudes presented in Eq. (5) 



the couplings G absorb the sign on imaginary z's as shown against the vertices in the 
figures mentioned above. The subscripts of G indicate the diagram and, if appropriate, 
the propagator gauge boson. On the other hand the superscripts {i or j) indicate whether 
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this coupling is arising at a vertex with an outgoing (i) or a propagator (j) gaugino. A 
superscript n in appropriate cases indicates the eigenstate of the scalar propagator. Cou- 
plings with superscript A and B are related by hermitian conjugation at the Lagrangian 
level, while involving the same set of fields in respective cases. There are two A-tjpe 
couplings for the t-channel processes in diagrams [Id), [Id') and le. Ai is associated 
with the vertex with the incoming e~ while A2 is associated with the vertex where the 
outgoing sfermion appears. Note that while A-type couplings appear in both s- and t- 
channel diagrams involving both final state gaugino and sfermion, 5-type couplings only 
appear in t-channel and associate with the incoming coupled either to a final state 
gaugino or one in the propagators. Also, note that in the most general case all the indices 
on the couplings G^'^ are nontrivial in the sense that couplings of a particular type with 
a generic structure may assume different values not only from diagram to diagram but 
also within a diagram and hence need special care while defining them. 

The qqg coupling follows from Eq. (C89) of Haber and Kane [0]. Special care should 
be taken in reading out the s-channel 7 and Z-mediated amplitudes of diagrams Figs. 2a 
and 2b when a single factor of e/sw should be replaced by the strong coupling constant gs- 
In addition one has to include in the amplitudes, the Gell-Mann matrices which appear 
in the squark-quark-gluino couplings. 

In all the couplings defined above we have taken the neutralino-mixing matrix in 
the B — basis, the neutralino-mixing matrix A^' in the 7 — Z basis, the two chargino 
mixing matrices U and V, to be real as is appropriate in an analysis that conserves CP. 
Following are the couplings in details. 

A-tjpe couplings in s-channel, where s = a,b,c and V stands for 7 and Z boson: 

G'^^(^> = C^{A^,),PL + A(^,)2PR)MnL + C^{A^,)3PL + A,(^,),PR)MnR (9) 

yl-type couplings in t-channel, where t = d, d', e: 

^t(L%]i>) = C^ji){Ai(^j)i(^^R^i^-)PL + Ai(^j-)2{L,R,D)PR) 

B~type couplings which appear only in t-channel: 

^lIr) = Cl(^R){Bi(^j)iL{R)PR + Bi(^j)2L{R)PL) 

Gf = C^iBaPR + B,2PL) (11) 

The generic structures of the terms on the right hand side of the above equations are 
defined below in reference to three final states broadly categorized as final states with 
neutralino, chargino and that with gluino. 
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For final states with neutralinos, one has: 

r-^A r-iA riA riB -i 

_ O2 — ^L{R) — ^L{R) — ~^ 

M{j)2L = Anj)lR = Bi(j)2L = Bi(j)iR = (12) 

For charged sleptons and d-type squarks (sneutrinos and w-type squarks) in the final state: 
Aim = V2lswQ,Kf^n + (^3,/ - Qfslr)Nly^,] 

Aim = Ao„ = -^r' with At") = ^ (^) 
^ ' ^' V2mw cos/3 Vsin/3y 

^,0-)4 = -V2(swQfNl^^),-sl,QfNl^^),) (13) 
For final states with charginos, one has: 

A{j)lL = Vilj)i , A{j)2L = A{j)lR = AU)2R = 

BjiiL,R) = Aji^L,R) with/ = 1,2 (14) 

For charged sleptons or d-type squarks (sneutrinos or w-type squarks) in the final state: 

ct = -cf = +i C| = C^ = C,^ = -C^^ = -1(+1) 
Cf = — —1 (appears only for the case of l'^ and d) 



A{j)i = Ui(j)i {Vi^j)i) , A. 



i{j)4. = u 

A,, - ^/ J with r<^rr"l - ( 

Ao)3 = with r,^(r,") = ^ f ^) (15) 

V2mvF cos/? ysin/?y 

Sj-; = with 1^1.. A 

Bii — Vii , — Bis — — (appears only for the case of l'^ and 0?) 

Finally, for gluino final states, we have: 

Al,^Al,^V2, A^^-A^^O (16) 

In our analysis, we will impose kinematic constraints prohibiting on-shell limits for 
the (sfermions and/or gauginos) propagators that would have led to resonances when the 
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sparticle widths, i.e. the Breit-Wigner form of the sparticle propagators, are not exphcitly 
used. This is quite justified as long as we are interested in studying situations reasonably 
away from such thresholds. But including sparticle widths in the propagators is a must 



when we are rather close to the thresholds; see for instance Ref. [^. We will demonstrate 
this effect separately in the study of pair-production of muons along with two neutralinos 
in a 4-body final state in e^e~ collisions with an emphasis on the role of finite widths of 
sparticles close to such thresholds. 



2.3 The differential cross sections 

The analytical expressions of the amplitudes squared of the various processes are quite 
lengthy and not very telling^. In the Appendix, we will simply give the amplitude squared 
of the sum of the three diagrams Fig. la-c, the contribution of which is the dominant 
one as will discussed later and explicitly shown in section 4. In principle, close to the 
\/s ~ 2mj threshold, diagrams (la) with an (almost) resonant intermediate sfermion 
should give the dominant contribution. Below this threshold, an important contribution 
will come from diagram (lb) with a virtual fermion, and in the case of associated pro- 
duction of sfermion-fermion-chargino final states, from diagram (Ic). However, these 
diagrams are not gauge invariant by themselves, and only the sum of the three contri- 
butions (la+lb+lc) is gauge invariant^. We will further use the simplification that the 
off-shell sfermion is of the same species as the real one, which is a very good approx- 
imation in the case of first and second generation sfermions where the mixing is very 
small and the Z-boson does not couple to different chiral sfermions. In the case of third 
generation sfermions, the mass splitting between the sfermion eigenstates is large if the 
mixing is strong and the contribution of the heavier sfermion is suppressed by its large 
virtuality. So here also, we carry on with the lighter species in the final state only. [Note 
that the case of first and second generation squarks is similar to the smuon case]. 



^The Fortran code which calculates the total cross section and includes all the amplitudes squared 
and the interferences can be obtained from datta@lpm.univ-montp2.fr. 

•^In a covariant gauge, the gauge dependent term is present only in the longitudinal component of the 
photon and Z-boson propagators and should therefore give a vanishing contribution when it is attached 
to the initial leptonic tensor in the limit of zero electron mass. However, in other (non-covariant) gauges 
such as the axial gauge, the gauge dependent terms are explicitly non-zero in the individual contributions 
of diagrams (la), (lb) and (Ic) and only the sum of the three contributions is free of this gauge dependent 
terms. In particular, diagram (Ic) which, for associated smuon and LSP production as will be discussed 
later, should give a contribution that is proportional to the small ©(Ml//!^) Z~x~X coupling in the 
gaugino region, can be neglected only when the calculation is performed in a covariant gauge. This is 
not automatically the case in a general gauge and to be consistent, we will include also this contribution 
in the analytical expression of the differential cross sections. We thank Ayres Freitas and Peter Zerwas 
for a discussion on this point. 
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3. The associated production cross sections 
3.1 The weight of the various contributions 

Before presenting our numerical results, let us discuss the contributions from diagrams 
Fig. la-e to the production cross sections, focusing on the two extreme cases where the 
higgsino mass parameter /i is much smaller or much larger than the gaugino mass param- 
eters M2. All through this analysis, we also presume the unification of gaugino masses at 
the GUT scale, leading to the relations Mi = ftan^ ewM2 ~ \M2 and m-g 3.5M2 
at the weak scale. 

If ^ Ml 2, the lightest chargino xf and the two lighter neutralinos x5 2 gaugino- 
like states with masses m^o ~ Mi, m^o ~ m^± ~ M2 ~ 2Mi, while heavier charginos 
and neutralinos are higgsino-like states with masses m^o ~ m^o ~ m^± ~ |//|. The 
couplings of the gauginos to electron-selectron or electron-sneutrino states are maximal, 
while the couplings of the higgsinos to sfermions and massless fermions vanish. [Note 
that in this case, the couplings of right-handed sfermions to the gauginos xf and X2 
also vanish]. In this case, only the contributions of the diagrams with light gauginos 
in the final state, and light gauginos exchanged in the diagrams Fig. Ic-le, have to be 
considered. [Contributions with heavier higgsino final states, in the case of smuon and 
sneutrino production, are suppressed, in addition to the smaller phase space, by one, two 
or even three powers of the small higgsino-lepton-slepton coupling.] In the case of r 
sleptons, the mixing can be very strong for large values of tan (3 and 11 leading to non- 
negligible higgsino-f-r or higgsino-f-z/,- couplings. In this case not only the s-channel 
diagrams (la) and (lb) would contribute but also diagram (Ic) for higgsino final states, 
although it might be suppressed by phase space. The situation for third generation t and 
h squarks is similar to that of the f slepton. 

In the opposite situation, |^| <^ Mi^2, the lighter chargino x^and neutralinos X12 "^^i^ 
be higgsino like and degenerate in mass. Their couplings to selectrons and sneutrinos will 
be very tiny and therefore the diagrams (Id, Id') and (le) will give very small contributions 
with light chargino and neutralinos in the final or intermediate states. In this case, 
processes involving heavier neutralino and chargino final [or intermediate in the case of 
Fig. Id-le] states may significantly contribute to the cross sections, but since these states 
are rather heavy, the processes will be kinematically disfavored [in the case of intermediate 
gaugino states, they will be suppressed by the larger virtuality of the particles too]. Again 
diagram (Ic) will give a very small contribution, but this time, because of the tiny final 
higgsino coupling to sfermion-fermion pair [except in the case of the top squarks because 
of the large value of m^, and f, b sfermions for very large tan /3 values]. 
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In the mixed gaugino-higgsino region, i.e. for |/i| ~ M2, and for associated production 
of third generation f leptons and t, b squarks with chargino and neutrahnos, the situation 
is more comphcated [as is well known] because of the large mixing in both the sfermion 
and gaugino sectors and all types of diagrams involving all gaugino states have to be taken 
into account. However, for a given final state, the cross section will be in general smaller 
than the one in the pure gaugino case, despite a more favorable phase space [since all 
neutralinos and charginos will have comparable masses in this case]. This is because the 
couplings of chargino and neutralinos to sfermions [in particular to sleptons] are reduced 
compared to the gaugino case. 

Thus, we will concentrate in our analysis on the more favorable case where the higgsino 
mass parameter n is much larger than Mi 2, i.e. on the gaugino LSP case, although we 
will give some illustrations in the other cases. 



3.2 Numerical analysis 

In this section we will illustrate the magnitude of the total production cross sections for the 
various processes discussed above. We will set for definiteness tan /3 = 30, and in the main 
part of the analysis, we will choose a large value of the higgsino mass parameter, fi = 500 
GeV, and a low value of the wino mass parameter, M2 ~ 100 GeV, leading to gaugino- 
like lighter chargino and neutralinos with approximate masses m^± ~ m^o ~ 2m^o ~ 100 
GeV, i.e. values close to the experimental LEP2 bounds [jT3|. The large tan/5 and /i 
values will imply a large mixing in the f and b sectors, leading to relatively light states 
compared to ft and q, respectively. For f sleptons, the no-mixing case is realized from 
the fi case, where the large value of tan f] and /i do not have a substantial effect since the 
smuon mass is too small. The situation for the two sneutrinos Ur and z>^ will be of 
course identical. 

For the soft-SUSY breaking scalar fermion masses at low energies, we use a universal 
value rrif =171?. To obtain the masses m? and one has to include the "D-terms" 
but since their contributions are rather small for of the order of a few hundred GeV, 

the mass eigenstates /i and /2 will be almost identical to the current eigenstates and mass 
degenerate. The degeneracy is lifted in the case of third generation charged sfermions by 
the mixing which is proportional to mf{Af — yutan/5) for the isospin — | b and f states 
and rrit^At — ji/ tan/3) for the isospin +| top squarks. In the analysis, we will set the 
trilinear couplings Af, which has no major effect in this context, to zero. 

Note that we have thoroughly cross checked our results against the corresponding ones 
from the program CompHEP |T6| in the cases of jlL,R-, i^fi and found very good agreement. 



Fig 3 shows the cross sections for the production of smuons and their sneutrino partners 
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in association with the hghtest neutrahno (3a) and the hghtest chargino (3c) as well as 
the production of the lightest fi slepton in association with the lightest chargino and 
neutralino (3b). The cross sections are shown for a cm. energy of 500 GeV as a function 
of slepton massQ, for the SUSY parameters discussed in the previous section. Shown are 
the summed up cross sections for the processes e'^e~ — > iix and the charge conjugate 
process e"'"e~ —>■ U*Xi which means that if the production of only one state is to be 
considered, the cross sections should be divided by a factor of two. 

In Fig. 3a, one sees that the cross section for the production of the LSP with the right- 
handed smuon fiR is larger than the one for the production with the left-handed smuon 
jiL and its partner sneutrino z/^, in particular near the kinematical threshold = 2m^. 
For instance, the cross section for the process e"'"e~ — > fifiRXi is close to 0.05 fb for 
m^^ = 260 GeV. This means that with the very high integrated luminosities, J C ^ 500 
fb~^, expected at future linear colliders [Q, 25 of such events can be collected within one 
year, opening the possibility of discovering for masses of 10 GeV above the threshold. 
The cross sections decrease quickly with larger m^, dropping to below 0.01 fb for all three 
processes [i.e. 5 events for the above luminosity] for masses close to 300 GeV. 

The production cross sections of the and its partner sneutrino in association with 
the lightest chargino are shown in Fig. 3b [fiR does not couple to the xf states in this 
regime], and as can be seen, they are approximately a factor of two larger than the 
corresponding production cross sections in association with the LSP, a mere consequence 
of the stronger charged current couplings compared to the neutral current ones. 

The cross sections for the associated production of the lightest f lepton with the 
lightest charginos and neutralinos. Fig. 3c, are similar to those of /il in spite of sfermion 
mixing. As mentioned earlier, the cross sections for the associated production of both 
types of sneutrinos, z/^ and Ur, with charginos is the same. 

The cross sections for the same associated production processes but for a center of 
mass energy ^/s = 1 TeV are shown in Fig. 4 as a function of the slepton masses for 
the same set of inputs as in Fig. 3. The production rates are smaller than in the pre- 
vious case, a consequence of the fact that the dominant contributing channels are the 
s-channel diagrams, which scale like 1/s at high-energy. This drop should however be 
partly compensated by the increase of the luminosity with energy [as expected for the 
TESLA machine for instance], allowing for a reasonable number of events for slepton 
masses a few tens of GeV above the kinematical threshold. 

'^In all the plots, we start the variation of the sfermion masses at approximately 2 (4) GeV above the 
kinematical threshold for y/s = 500 GeV (1 TeV). These values typically correspond to the total decay 
width of the produced sfermion with masses of 250 (500) GeV. Close to the kinematical threshold, a more 
dedicated analysis is required as will be discussed in section 4. 
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Figure 3: Associated production cross sections for sleptons and gauginos in e+e" collisions 
for a cm. energy s/s = 500 GeV as a function of the slepton masses. For the SUSY 
parameters, we have set m^^ — rn^^, M2 — 100 GeV, = 500 GeV and tan /3 = 30. 
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Figure 4: Associated production cross sections for sleptons and gauginos in e+e" collisions 
for a cm. energy ^/s = 1 TeV as a function of the slepton masses. For the SUSY 
parameters, we have set m^^ — rn^^, M2 — 100 GeV, = 500 GeV and tan /3 = 30. 
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Fig. 5 shows the associated production cross sections in the smuon sector in two 
particular cases for a cm. energy ^/s — 500 GeV. In Fig.Sa, the cross sections for p,L and 
i>i^ production in association with the heavier neutral gaugino X2 displayed [again, 
jiR does not couple to X2 this case]. Approximately, they are the same and equal 
to the ones with LSP final states close to the kinematical thresholds, despite the fact 
that m^o ~ 2m^o. In Fig. 5b, we show the cross section for a fixed slepton mass of 
mi = 275 GeV and a varying M2 parameter. While the variations of a{e~^e~ — > ^jliXi) 
and a{e^e~ I^'P'LxV) mild, the cross section a{e^e~ ^fiP'Lxt) drops quickly with 
increasing M2. Therefore, the cross sections for the associated production of smuons with 
heavier chargino states will be in general smaller than what was shown in Figs. 3 and 4. 
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Figure 5: Associated production cross sections for sleptons and gauginos in e'^e" collisions 
for a cm. energy ^Js — 500 GeV as a function of the slepton masses with x = X2 (^) 
as a function of the gaugino mass parameter M2 with X — X12 (b) ■ 
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Finally, Fig. 6 shows the production cross sections for the case where the hghtest 
neutralino is higgsino-hke in Fig. 6a [with fj, — 100 GeV and M2 — 500 GeV] or a mixture 
of higgsino and gaugino states in Fig. 6b [with M2 — ji — 150 GeV leading to xt masses 
m^± ~ 2m^± ~ 200 GeV]. As mentioned earlier, the cross sections are extremely small 
in the higgsino case and the situation should be the same for the production of smuons 
with X2 , and also for the production of f's with higgsinos for not too large values 

of IX and tan/?. However, for the mixed case, the cross sections are only slightly smaller 
than in the gaugino case. This is due to to the fact that all neutralino and chargino states 
[which share the gaugino couplings] have to be taken into account, leading to an overall 
contribution which is not very far from the one for the pure gaugino case. 
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Figure 6: Associated production cross sections for sleptons and gauginos in e^e" collisions 
for a cm. energy ^ — 500 GeV as a function of the slepton masses for the higgsino case 
(a) and the mixed gaugino-higgsino case (b). 
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Let us now turn to the case of t and b squarks. The cross sections for the associated 
production of these squarks with the hghtest neutrahno and chargino are shown in Fig. 7 
for center of mass energies of 500 GeV (a), 1 TeV (b) and 3 TeV (c). In the 500 GeV 
case where phase space plays an important role due to the heaviness of the top quark, 
the largest cross section is obtained in the case of e~^e~ bbiXi but it barely reaches 
the level of 0.03 fb for m^^ ~ 260 GeV, leading to not much more than 10 events for 
a luminosity of / £ ~ 500 fb^^. It is followed by the cross section for e~^e~ — > btixf 
which is approximately a factor of two smaller. The cross section for e"^e~ ttiXi is one 
order of magnitude smaller because of the strong suppression by the phase space [since 
for m^^ > 250 GeV, rrit + m^_^ + m^o > 475 GeV, i.e. close to y/s = 500 GeV] while the 
process e~^e^ tbixf is not kinematically accessible at this energy. 

At higher energies. Figs. 7b and 7c, the largest cross section is a{e^e~ biixf), which 
at y/s = 1 TeV can reach the level of 0.1 fb for a stop mass of a few tens of GeV above 
the kinematical threshold, followed respectively, by the cross sections a{e~^e~ tiiXi), 
a{e~^e~ — > tbixf) and a{e~^e~ — > bbiXi) which, up to a factor of two, have the same size. 
This hierarchy of cross sections near the threshold follows from the dominance of final 
states with stop squarks since the main contribution comes from the s-channel production 
of stops with photon-exchange [because of the larger Q^] and by the fact that charged 
current processes are more favored. At even higher energies. Fig. 7c, the trend is similar 
to the previous case, but the cross sections are smaller since the contributions of the 
dominant s-channels scale like 1/s. 

Finally, the cross sections for the associated production of third generation squarks 
with gluinos are shown in Fig. 8. Here, there are two kinematical situations which are 
relevant for a (genuine) three-body particle final state: (i) if the gluino is heavier than 
squarks and in this case, it is the only way to produce gluinos in e~^e~ collisions [except for 
the loop induced production of gluino pairs [0 for which the cross sections are too small] 
and (ii) if the gluinos are lighter than squarks, with the latter not being kinematically 
accessible in pairs in e~^e~ collisions [otherwise, one can produce squarks which then decay 
into a quark and a gluino pair, with a branching ratio that is dominant since it is a strong 
interaction process]. The cross sections for gluino production in association of bi and ti 
squarks are shown for cm. energies of 1 TeV (Figs. 8a, 8c) and 3 TeV (Figs. 8b, 8d) in 
these two kinematical regimes. In the first situation [gluino heavier than squarks], the 
production cross sections can reach the level of a fraction of a femtobarn leading to a 
few tens of events for the expected luminosities, / £ ~ 500 fb~^. In the second situation 
[squarks heavier than gluinos], the cross sections are rather tiny [below 0.01 fb] and higher 
luminosities are needed to detect these final states. 
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Figure 7: Associated production cross sections for squarks and electroweak gauginos in 
e+e" collisions for a cm. energies of ^/s — 500 GeV (a), 1 TeV (b) and 3 TeV (c) as a 
function of the squark masses. For the SUSY parameters, we have set m^^ = ''^qn-i ^2 — 
100 GeV (150 GeV in c), = 500 GeV and tan/3 = 30. 



21 



0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 h 



' 1 
(a) 


(7(e 


1 ' 

e qqg) [Vo\ - 




V~s 


= 1000 GeV 


\ 


9 


= 350 GeV 


\^ \. 

N 

N 








N. 




1 







0.025 



0.015 - 



0.005 - 



0.4 
0.35 
0.3 
0.25 
0.2 
0.15 
0.1 

H 0.05 h 




350 



300 



0.05 



0.04 



0.03 - 



— I ' 1 ' 1 

(b) cr(e+e- ^ qq~g) [fb] 

= 3000 GeV 
m-n = 700 GeV 




hh^g 



400 



500 



600 



700 



— I ' 1 ' 1 ' r 

(d) cr(e+e- ^qqg) [fb] 

= 3000 GeV 
m, = 525 GeV 




0.02 - \tUg 



0.01 K ~ 

\bh~g 



540 560 580 600 



1500 1600 1700 1800 1900 2000 



Figure 8: Associated production cross sections for bi and ii squarks and gluinos in e'^e" 
collisions for a cm. energy ^/s — 1 TeV (a,c) and 3 TeV (b,d) as a function of the squark 
masses. For the SUSY parameters, we have set M2 = 350 GeV (525 GeV) in a-c (b-d), 
H = 500 GeV and tan (3 = 30. The increase of the cross section in (a) for sbottoms for 
^fti ^ '^a ^'^ small virtuality of bi in this case. The kink in (b) is due to the 

heavier 62 (exchanged in the s-channel digram with a Z) becoming close in mass to the 
gluino. For stops, this is not apparent in the curves because of the large value of m^. 
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4. The cross sections in various approximations 

As mentioned in section 2.3 and 3.1, the bulk of the production cross sections for the 
associated production of sfermions, fermions and electroweak gauginos is due to the con- 
tributions of the s-channel diagrams, Figs, la-lc. In particular, near the kinematical 
threshold for sfermion pair production, ^/s ^ 2171 j, the contribution of digram (la) with 
an almost resonant sfermion is dominant. However, this diagram is not gauge invariant by 
itself, although in a covariant gauge, the gauge parameter dependence drops out. Further- 
more, very close to the kinematical threshold, the zero-width approximation is not valid 
anymore and one has to include the total decay width of the sfermion in the Breit-Wigner 
form of the resonant propagator. This is performed by substituting in the sfermion prop- 
agator, the squared mass by m^ — im^Fj. In fact, near this kinematical threshold, the 
three-body approximation is not accurate anymore and for a more reliable result one has 
to consider the four-body process with the production of two off-shell sfermions which 
then split [or decay if one is above threshold] into fermions and gauginos. Here again, the 
sole contribution of this double resonant diagram [which, in general, is the only one taken 
into account in experimental analyses, see Ref. |]18| for instance] is not gauge invariant and 
one should take into account additional diagrams with single resonant states [i.e. with 
only one intermediate sfermion state or with neutralino or chargino states] ||11| |. 

In this section, we will therefore discuss in more details the magnitude of the pro- 
duction cross sections in the various approximations taking as example the associated 
production of the right-handed smuon, a muon and the LSP, e+e~ Ai^/i_RXi cind the 
charge conjugate final state ii~jl\xi- The total decay width of the right-handed smuon, 
assuming that the LSP is a pure bino [which implies that jiR does not couple not only to 
the heavier charginos and neutralinos which are higgsino-like but also to the gauginos xf 
and is identical to the partial decay width into a muon and the LSP which is given 
by [the couplings given in section 2.2 have been simplified] 

= \4K. - "-^N[Xm,, (l - (17) 

For a 250 GeV fin, assuming as usual M2 = 100 GeV with tan /3 = 30 and ^ = 500 GeV, 
the total decay width is F^^ ~ 1.15 GeV. 

In Fig. 9a, we show the cross section for this process as a function of the cm. energy 
[varied from 400 to 500 GeV] for ^ = 250 GeV and the SUSY parameters given above, 
retaining F^^. While the short-dashed curve shows the contribution of the diagram (la), 
the long-dashed curve shows the contribution of the three s-channel diagrams (la), (lb) 
and (Ic), and the full line shows the total contribution including all diagrams. 
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Figure 9: The total cross section for the associated production process e~^e~ fJ^'finXi 
as a function of the center mass energy in the exact case (sohd hues) and in the approx- 
imations of the gauge invariant set of diagrams (dashed hues) and the resonant smuon 
diagram (dotted hues). In panels (a) and (c) the total smuon width is included, while in 
panel (b) it is set to zero. 
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As one can see, it is a very good approximation [in the present case] to include only 
the gauge invariant set of contributions (diagrams la+lb+lc), since the full line and the 
long-dashed line are almost overlapping in the entire range of ^/s values shown in the 
figure. Far below the kinematical threshold, the contribution of diagram (la) is rather 
small, but closer to the threshold where the exchanged smuon is almost on mass-shell, 
it becomes the dominant one. This is shown in a more explicit way in Fig. 9b, where 
we zoom on the threshold region. A few GeV near y/s = 2m^^, the contribution of 
the almost-resonant smuon exchange diagram represents 95% of the total cross section. 
In Fig. 9c, we display the effect including the total width of the smuon in the resonant 
diagram (la). As can be seen, the production cross section is slightly lower [by ~ 5%] 
compared to the zero-width case. Fig. 9b. Therefore, close to the kinematical threshold, 
one has to include all the s-channel diagrams as well as the finite width of the smuon to 
have a rather accurate estimate of total cross section. 

Let us now turn to the discussion of the four-body production process, e~^e^ 
/i^xV Xi- Taking into account only the gauge invariant set of the double-resonant 
diagram, e+e^ 1^*rP'*r ~^ plus the single resonant ones with either smuons 

or gauginos in the intermediate states, 10 Feynman diagrams [the generic form of which 
can be obtained from the set of diagrams (la), (lb) and (Ic) of Fig. 1, where the smuon 
splits into a muon and a neutralino] are to be considered. We have calculated the total 
production cross section, including the finite widths of the smuons, in this approximation 
using the package CompHEP [|T^ . This should be a good approximation^ corresponding to 
what occurs in the three-body production process. We have compared our results with 
the ones obtained in Ref. |TT| and found perfect agreement. 

The four-body cross section, as a function of the cm. energy, is shown by the solid 
line of Fig. 10 for the set of inputs chosen previously. It is compared with the cross section 
for the three-body process below and above threshold^ [long-dashed lines] and with the 
cross section of the two-body e^e~ i^^i^li production above the threshold with the 
subsequent decays of the smuons into /xxi final states [short-dashed lines]. The on-shell 

^The full calculation of the four-body final state process involves 84 Feynman diagrams. Even if one 
neglects the diagrams involving the couplings Z-x^-yp^ ^-^i?-X2 [which vanish in the gaugino LSP limit 
in a covariant gauge] and Z-fli-fl2 [which is zero in the no-mixing case] , 44 Feynman diagrams remain. 
Even with CompHEP, the calculation in this simpler case was consuming an enormous amount of time and 
we did not perform it. This is partly due to the huge volume of the kernel matrix element squared and 
partly to the Monte-Carlo integration on the 4-body final state which calls for a large number of points 
and iterations to have a better convergence of the integral. 

Note that for the three-body process, we include for the below threshold usual, both the 

cross section for fiji and for the conjugate state /ijj. Above the threshold, we have to divide the cross 
section by a factor i to take into account the fact that, after the decay of the (on-shell) smuons, we have 
two identical neutralinos in the final state. 
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smuons are assumed to decay 100% of the time into fiXi final states. As can be seen, below 
the kinematical threshold the four-body and the three-body cross section are almost 
overlapping, which means that it is really a very good approximation to consider in this 
regime the much simpler three-body production case. Sufficiently above the kinematical 
threshold, the three-curves also overlap as it should be. A few GeV above the threshold, 
the three-body and two-body approximations fail badly, and one has to consider the full 
process with two virtual smuons including their total decay widths in the propagators]]. 

In fact, for a better description of the cross section around the kinematical threshold, 
which is needed to meet the experimental possibility of measuring the produced smuon 
mass with a scan near threshold with a precision better than 100 MeV 0, one has not 
only to use the full four-body production process, e+e~ including the 

smuon finite widths, but also to take into account some refinements such as Coulomb 
re-scattering effects, initial state radiation from the emission of collinear and soft photons 



and beamstrahlung effects, as discussed in Ref. |TT 
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Figure 10: Cross section for the two-, three- and four-body processes leading to the final 
state with intermediate right-handed smuons, around a cm. energy of 500 

GeV. For the SUSY parameters, we have set ~ 2Mi = 100 GeV, /i = 500 GeV and 
tan/3 = 30. For the three-body curve, the kink at y/s ^ 2(mj + r^) is due to the "brutal" 
inclusion of the factor | to avoid double counting above threshold; the approximation 
slightly above threshold is thus rather bad (see text). 



^Of course, one has to consider the fuU contribution of the gauge invariant set discussed above. For 
•^i = 485 GeV, this gives a total cross section of 0.046 fb for the set of inputs of Fig. 10, while the 
inclusion of only the doubly resonant diagram gives a cross section of only 0.037 fb. 
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5. Conclusion 



In this paper, we have analyzed the production at future e^e~ colhders of second and 
third generation sleptons and third generation squarks, in association with their partner 
leptons and quarks and charginos or neutrahnos [and gluinos in the case of squarks], 
e^e~ —>■ f fXi in the framework of the MSSM. Analytical expressions for the differential 
cross sections have been given in the approximation where only the s-channel contributing 
diagrams are taken into account, an approximation which has been shown, a posteriori, to 
be very good. Taking the example of right-handed smuon production, we have shown that 
consideration of only the three-body processes is a very good approximation below the 
kinematical two-sfermion production threshold, but that close to these threshold the full 
four-body process e^e~ ffxXi including the finite total decay width of the sfermions, 
have to be considered |Tl[ . 



We have shown that some of these three-body processes can have sizeable enough 
production cross sections to allow for the possibility of discovering sfermions with masses 
a few tens of GeV above the kinematical threshold for sfermion pair production ^/s = 2m j 
in favorable regions of the SUSY parameter space. This is due to the very high luminosities 
expected at future linear e~^e^ colliders, J C 500 fb~^, i.e. three orders of magnitude 
larger than in the case of LEP2 where the experimental limits on the sfermions masses 
do not exceed the beam energy. Associated production of squarks and gluinos offers, for 
some range of sparticle masses, the unique direct access to the gluinos in e+e" collisions, 
if they are (slightly) heavier than squarks. 

The final states discussed in this paper should be clear enough to be easily detected 
in the clean environment of e~^e~ colliders. However, for a more precise description, 
detailed analyses [which are beyond the scope of the present paper] which take into account 



Standard Model and SUSY backgrounds [as those performed in Ref. [jTTI] in the case of 
smuon production near threshold] as well as higher order effects and detection efficiencies, 
have to be performed. 
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Appendix: the differential cross section 

In this appendix, we present the analytical expression of the differential cross section for 
the associated production of sfermions with charginos or neutralinos 



e+(pi) e (P2) fh{P3) XiiPd /'(Ps) 



(A.l) 



where h — L,R\s the handedness of the produced sfermion [later on, we will use h' — R,L 
which is the other possible helicity]. To simplify the expressions, we will work in the 
approximation in which the accompanying fermion in the final state is massless, leading 
to zero-mixing in the sfermion sector [these expressions are therefore not accurate in the 
case of associated production with final state top quarks]. Furthermore, we will take 
into account only the contributions of the gauge invariant set of diagrams Fig la, lb and 
Ic with s-channel exchange of photons and Z bosons which, as shown in the numerical 
analysis in section 4, is a very good approximation to the full cross section. 

The complete spin-averaged matrix element squared for vanishing final state fermion 
mass and zero-mixing in the scalar sector, is given by 

(47rQ;)3 



AN, 



c 2 



W 



i>j=a,b,c 



(A.2) 



where Nc is the color factor and a is the fine structure constant. The Tij are the amplitudes 
squared with i = j = a,b,c of the contributions of diagrams la, lb or Ic where the 7 
and Z boson amplitudes have been added, and the interference terms arc obtained for 
i 7^ j = a,b,c; here the effects of permutations of suffixes are already included in the 
original terms and hence such permutations are to be left out. 
The diagonal terms are given by: 



T = 

n.n. 



42 



0} ^ 

{s - m\Y s^c^s(s - m|) 



zff 



Tr 



bb 



42 



{{q - p,y - m'jy 



Q) , {al + vl){af-Vff , 2Qf{af-Vfy 



. s2 ^ Sy^Cy^{s - m|)2 



+ 



2ae{af-Vf) j Ve{af-Vf) 



.^w'^wi^ ^z) 




(A.4) 
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Ql 



E 



42 



2 2 

m^. — m? 

Jh 



)K^ + K'^ + K'i} 



28 



2 2/ 2 



while the non-diagonal terms are given by: 

Tab 



42 



+ 



^ 



2 2/" 2 



zff) \ ^'12 



K" 



(A.6) 
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QxiQi 

s I s 



13 



^w'^wi^ ^z) 



— --2 2 



(A.7) 



1 



6c 



E 



42 



sis 



23 



+ ^23) + 
+ 



1 



s^ciy[s-mi) sj ^ s J 



The factors K, which involve the dependence on the momenta, read: 





= P45 (4^13^23 -"^/^«) 


K2 




K 


= P46(Pl5 -P25)S 




= P14P25 +P1BP24 




= (P13P24 +Pl4P23)P3B 




= TT^Xi^XiSPs6 
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K2 = S {{p,, - p^,)Ps^ - ip,^ - ^24)^35} 

= K12 - 2p,,p^^ S 

Kl3 = 2 {{p,,p.2, + ^15^23)^34 - (P13P24 + ^14^23)^35} - (4P13P23 - TTT'X S)P45 

K3 = "^Xi"^X. {2(Pl3P25 +P15P23) -P35S} 

^23 = 2(p,3P2B +PibP23)P45 

^23 = "^Xi^^X, {«(Pl6 +P25) -4P1BP25} 

^23 = {« (Pi4 + P24 ) + 2 (P14P2B + P16P24 ) } P36 - 4P1BP25P34 

^23 = ^23-2sp3BP4g 

^23 = '^P2,{Pl3P4, - Pl4Ps,) + + 

Kg = sm^m^^{p,,-p,,) 

= s{p,,-pjp,, (A.9) 

In the above expressions, h is the suffix for the chirality (L or R) of the final state sfermion 
and h' is the other chirality {R or L); the latter suffix appears only on the couplings O^. 
q — Pi +P2 is the sum of incoming 4- momenta. We used the notation 2ve = C£ + and 
2ae = C*! — and similarly for and Uf. The correspond to the scalar products 
Pi -pj and s — (f — (pi +^2)^- All other couplings and parameters are defined in the text. 



Note that the antisymmetric part in each of the above terms, proportional to ae are 
pulled out apart and when, integrated over the full angular domain the corresponding 
contribution which antisymmetric, vanishes. 

The differential cross section is obtained by dividing by the flux and multiplying by 
the phase space. 

The integral over the phase space, to obtain the total production cross section, is then 
performed numerically. 
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